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1 .   Introduction 

Hoskins  and  Bretherton  (1972)  have  shown  analytically,  and  Williams 
(1967,  1972)  has  shown  numerically  that  discontinuous  fronts  can  form  within 
a  finite  period  of  time  if  no  turbulent  diffusion  is  present.   These 
studies  suggest  that  a  discontinuity  will  form  within  24  to  36  hours  when 
reasonable  initial  conditions  are  used.   If  turbulent  diffusion  is  present, 
it  can  be  expected  that  a  balance  will  be  achieved  between  the  fronto- 
genetic  advections  and  the  turbulent  diffusions  of  heat  and  momentum. 
The  front  should  remain  in  this  state  of  quasi-balance  as  long  as  the 
large  scale  deformation  field  causes  frontogenetic  advections  around  the 
front.   The  purpose  of  this  study  is  to  obtain  and  examine  steady-state 
frontal  solutions.  Welander  (1963)  obtained  steady-state  frontal  solutions 
with  a  boundary  layer  technique,  but  his  formulation  did  not  involve 
a  frontogenetic  large-scale  deformation  field. 

The  turbulent  diffusion  coefficients  for  heat  and  momentum  are  taken 
to  be  constant  in  this  study,  although  these  coefficients  are  actually 
functions  of  the  wind  and  temperature  fields.   Smagorinsky  (1963)  and 
Leith  (1968)  have  proposed  forms  for  the  horizontal  diffusion  coefficients 
which  are  dependent  on  the  local  wind  field.  Many  studies  have  treated  the 
form  of  the  vertical  diffusion  coefficients  near  the  surface  (Estoque  1960, 
Deardorff  1970,  1972).   In  this  study  we  will  carry  out  a  large  number  of 
integrations  with  different  values  of  the  various  constant  coefficients. 
The  range  of  values  will  include  the  mean  values  which  are  given  by  the 
various  nonlinear  formulations. 

This  investigation  will  employ  the  numerical  model  of  Williams  (1972) ; 
this  reference  will  hereafter  be  referred  to  as  W72,  The  physical  model  is 
essentially  the  same  as  one  of  the  models  treated  by  Hoskins  (1971)  and 
Hoskins  and  Bretherton  (1972).   The  frontogenesis  is  forced  by  a  non- 
divergent  horizontal  wind  field  which  contains  stretching  deformation. 
The  hydrostatic  primitive  equations  are  used  with  the  Boussinesq  approximation. 
In  the  model  the  time  dependent  quantities  are  functions  of  y  and  z  only. 
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The  basic  deformation  wind  field  is  constant  in  time  and  is  independent 
of  height,  except  in  the  surface  boundary  layer. 

In  Section  2,  the  basic  forecast  equations  are  developed,  and  a 
simplification  which  keeps  the  problem  two-dimensional  is  discussed. 
The  initial  conditions  are  presented  in  Section  3.   The  solutions  which 
contain  only  horizontal  diffusion  are  discussed  in  Section  4,  and  those 
which  contain  only  vertical  diffusion  are  discussed  in  Section  5.   The  com- 
bined experiments  which  include  both  types  of  diffusion  are  examined  in 
Section  6  and  the  conclusions  are  given  in  Section  7. 
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2 .   Basic  equations  ^    :  ,     .  m 

The  hydrostatic  Boussinesq  equation  with  diffusion  may  be  written: 

1^^  +  V  .  OK)  +  |-  (wV)  +  V(|)  +  ffe  X  V  =  AV2  V  +  ^m_|   ^     ^^.l) 

^   +  V  .  (ev)  +  ^  (we)  =  Av^e  +  c^  ^^ 


V  •  ^  +  ^  =  0,  (2.3) 


^   =    ge 


(2.4) 


where  the  notation  is  the  same  as  used  in  W72.   The  qxoantity  A  is  the 
horizontal  turbulent  diffusion  coefficient;  note  that  the  same  value  is 

used  for  both  heat  and  momentum.   The  quantities  C  and  C^  are  the  vertical 

m      6 

turbulent  diffusion  coefficients  for  momentum  and  heat  respectively. 
A  convective  adjustment  process  is  represented  by  the  function  Q  .   The 
boundary  conditions  at  the  top  and  bottom  of  the  domain  are: 

w(x,y,o,t)   =  w(x,y,H,t)   =  0 

V(x,y,o,t)  =  0 


C 


m 


^  (x,y,H,t)  =0  (    '  ^^'^^ 


B: 


z 


Cg^  (x,y,o,t)  =  Cg  ^  (x,y,H,t)  =  0 

^z  hz 

The  lower  boundary  conditions  are  realistic  for  fronts  over  land,  but 
the  upper  lid  boundary  conditions  give  only  a  rough  approximation  to  the 
tropopause.   Upper  frontogenesis  has  been  treated  by  Hoskins  and  Bretherton 
(1972)  and  Hoskins  (1972).   In  this  paper  our  main  concern  is  with  fronts 
in  the  lower  troposphere. 
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The  following  fields  which  satisfy  the  boundary  conditions  (2.5) 
are  approximate  solutions  to  the  steady-state  equations   [(2.1)  -  (2.4)] 

/        -Bz  ~Bz   ,     ,  . 

V  =  U  =  D  I  [x(l  -e    cos  Bz)  +  ye    sin  Bz]^ 

-[y(l  -  e    cos  Bz)-xe    sin  Bz]J  [ 

(2.6) 
w  =  0 

^  =  I  H  -D^(x^  +  y^)/2  -  fDxy 

e  =  0 

where  B  =  (f/2Cin)'^  and  where  D  and  f  are  constant.   These  relations 
satisfy  the  steady-state  equations  if  some  small  advection  terms  in  the 
boundary  layer  are  neglected.   These  small  terms  are  usually  dropped  in  the 
Ekman  theory  and  they  are  of  the  order  of  the  Rossby  Number  which  for 
these  fields  is  D/f.   In  this  paper  D/f  =  0.1  for  all  cases. 
We  now  subdivide  our  dependent  variables  as  follows: 

V  =  U  (x,y,z)  +  u(y,  e,  t)^  +  v  (y,z,t)  j^ 
w  =  w(y,  z,  t) 

e  =  e(y,  ^,  t)  ^    '  ^^'^^ 

<|>  =  f  (x,y)  +  7T(y,  z,  t) 

It  is  noted  that  all  departures  from  the  steady  deformation  solutions  (2.6) 
are  assumed  to  be  independent  of  x.   If  we  substitute  the  relations  (2.7) 
in  the  i,   component  of  (2.1)  we  obtain: 


^  +  ^  (uv)  +  S_  (wu)  +  uD(l  -  e"^^  cos  Bz) 

St    by  hz 

-D[y(l  -  e"^^  cos  Bz)  -  xe"^^  sin  Bz]  ^  +  vD  e"^^  sin  Bz  (2.8) 

-HwDS  [xe  ^(1-  cos  Bz)  +  ye"^^  sin  Bz]  =  fv  +  A^^u  +  CmS^u 

Please  note  that  the  assumption  that  u  is  independent  of  x  for  all  time 
is  violated  in  (2.8)  in  the  two  terms  where  x  appears.   This  indicates 
that  the  other  quantities  which  were  assumed  independent  of  x  are,  in  fact, 
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x-dependent  and  this  would  bring  other  terms  into  the  eqiiations.   In  this 
study  we  will  neglect  these  added  terms  and  apply  (2.8)  and  the  other 
equations  at  x  =  0,   If  the  initial  u,  v,  w,  0,  and  tt  fields  are  independent 
of  X  the  error  from  this  approximation  will  grow  slowly  and  it  will  be 
confined  to  the  boundary  layer  since  the  x-dependent  tetms  in  (2.8)  are 
zero  outside  the  boundary  layer.   Since  frontogenesis  occurs  very  rapidly, 
it  is  expected  that  the  development  of  x-variations  in  the  dependent 
variables  would  have  only  a  small  effect  on  the  resulting  quasi-steady 
front.   In  any  case  these  effects  could  not  be  observed  in  the  atmosphere, 
because  atmospheric  fronts  always  have  some  variation  in  the  basic  fields 
along  the  front  which  would  be  much  more  important. 
When  Eq.  (2.8)  is  evaluated  at  x  =  0  it  becomes: 

-Bz 
hu  +  h(uv)  +  _^wu)  -  Tu  ^-h  rv^  -  v^  e    sin  Bz 

St   ^      hz  ^     ^    ^ 

-VwS  (e  ^  sin  Bz)=fv  +  AS  u  +  GmS  u   , 

(2.9) 


Sz  ^2      ^^2 


-Bz 


where:   V  =  -Dy  and  p  =  1  "  e'  cos  Bz.  The  j  component  of  (2.1)  applied 
at  X  =  0  takes  the  form: 

-Bz 
bv  +  _^(w)  +  ^  (wv)  +  fS_(vV)  -  uSv  e    sin  Bz 

5t   Sy       Sz        ^        ^ 

-wva  (e'^^ces  Bz  )=  -^  -  fu  +  k^     +   Cj^v  ^  (2.io) 

^  ay        ay2      ^i2 

In  a  similar  manner  Eq.  (2.2)  may  be  written: 

aa  +  ^(v9)   +  ^(we)  +  rvSa  =  a^  +  cgS^e  +  q^  ,  (2. id 

St  ^  Sz  Sy  Sy2  Sz2 

and  Eqs.  (2.3)  and  (2.4)  become: 

Sv  +  Sw  =  0  ,  (2.12) 

^     Sz 

Sn  =  sa  .  (2.13) 

dz  BS 
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The  boundary  conditions  at  the  top  and  the  bottom  of  the  domain  are 
given  by: 

u  =  V  =  w  =  Cq_^  =  0,  z  =  0 

Sz  ,  (2.14) 

Cm^  =  Cm^  =  w  =  C0_^  =0,  z  =  H 
bz  bz  bz 

where  H  is  the  distance  between  the  rigid  plates.  The  no  slip  condition  at 

the  lower  boundary  and  the  no  strex  condition  at  the  top  roughly  represent 

the  troposphere. 

If  we  define  the  vertical  average  of  a  quantity  as 
H 


■>-i  l 


<()>=!   \  (   )dz, 
H   ^ 

and  integrate  the  hydrostatic  equation  (2.13)  with  respect  to  z  and  remove 
the  vertical  mean,  we  obtain 

z        z 

TT  -  <TT>  =  E-  \    edz   -  <  \edz  >     .  (2.15) 

eo    L  ^  o        -I 

Take  the  vertical  average  of  (2.12)  and  use  the  boundary  conditions  (2.14), 
which  gives 

^  <v>  =  0. 
by 

This  equation  states  that  the  total  mass  flux  in  the  j  direction  is  inde- 
pendent of  y.   In  W72  a  symmetry  argument  was  used  to  show  that  this  flux 
must  be  zero.   This  argument  does  not  hold  strictly  in  our  case,  but  it 
can  be  expected  that  at  large  distances  from  the  frontal  zone  the  distur- 
bance mass  flux  will  vanish. 
Thus  we  set 

<v>  =  0.  (2.16) 
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If  we  utilize  the  development  in  W72  we  can  rewrite  (2.10)  in  the  following 
form: 

hv     +  ^  [w  -  <w>]  +  ^(wv)  +  _^  (pvV  -  <:rv>V) 
Bt    ^  Sz      ^ 

-s„  r       ~Bz    ,      „  -Bz         .      „      ,  r  -s    ,   -Bz 

-   6V  lue        sin  Bz   -<ue  sm  Bz>J    -  Vjwo   (e  cos  Bz) 

Sy  .  L  Sz 

-<w^_(e"^^eoe  Bz)>l   =  -   h  (tj  -  <tt>)    -    f  (u   -  <:u>) 
hz  ■J  "^ 

+  kh^v  +  C^\  ^^v     +  I  /^\  1 

b^         [bz^       n[hz)^  ^  oJ  •  ^^•'^'^^ 

The  Eqs.  (2.9),  (2.11),  (2.12),  (2.15),  and  (2.17)  forma  complete 
set  which  can  be  solved  by  a  pure  marching  process.  The  finite  difference 
equations  conserve  mean  squares  in  the  advection  terms  and  are  described 
by  Williams  (1967).   The  frictional  terms  are  evaluated  at  the  previous 
time  step  to  avoid  computational  instability.   Computational  boundaries 
are  introduced  at  y  =  ±  Y  and  these  are  the  same  as  used  in  W72.  On  these 
boundaries  u  and  0  are  held  fixed  and  v  =  0, 

3.   Initial  conditions 

The  initial  conditions  are  the  same  as  those  employed  in  W72.   The 
initial  temperature  field  is  given  by 

e(y,  z,  o)  =  ^T  (z  -  h/2)  -  a(2/rr)  arctan  (sinh  ay),  (3.1) 

hz 

where  a  =  f  n  H   (gg^   ^Qi/hz)'^^ .      The  quantity  hQj/hz,   which  is  constant, 
is  the  initial  static  stability  and  a  is  one  half  the  total  horizontal 
temperature  variation 

The  initial  x  -  component  of  the  velocity  is  given  by 
u(y,  z,  o)  =  2  .gaa  (z  -  H/2)  sech  ay.  (3.2) 

TT   fGo 

This  equation  is  obtained  by  substituting  (3.1)  into  the  thermal  wind 
equation  and  integrating.   In  this  study  the  v  field  is  the  divergent  part 
part  of  the  wind.  We  initialize  the  v  field  by  equating  it  to  the  diver- 
gent part  of  the  wind  which  is  obtained  from  the  quasi-geostrophic  equations. 
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The  details  are  given  in  W72. 

It  is  noted  that  these  initial  conditions  do  not  satisfy  the  boundary 
conditions  on  u,  v,  and  0  at  the  upper  and  lower  boundaries.  A  period  of 
adjustment  will  be  required  to  form  the  surface  friction  layer  and  some 
oscillations  may  be  observed  later. 

4,   Numerical  experiments  with  horizontal  diffusion 

All  numerical  results  to  be  shown  use  the  following  values  for  the 
constants: 

f  =  lO"'^  sec"-"-  ,  £  =  .0327  m  sec-2(OK)-l 

Qo 

H  =  9  km,   Y  =  1800  km  >  '  (4.1) 

D  =  lO'^sec""""  ,   oej/^2  =  ^°^   ^'■^ 
a  =  12.56K 

These  values  were  also  used  in  W72. 

In  this  section  we  examine  the  numerical  solutions  which  include  only 
horizontal  diffusion,   so  that  Cm  =  Cq  =  0.  Many  numerical  experiments 
are  examined  in  this  section  which  have  values  of  the  horizontal  diffusion 
coefficient  ranging  from  A  =  10""  m2sec~-'-  to  A  =  lO^m^sec"-'-,  All  experiments 
in  this  section  have  a  vertical  increment  of  Az  =  333  m.   The  horizontal 
mesh  length  is  either  20  km  or  60  km  depending  on  the  expected  frontal 
scale.   The  convective  adjustment  process  Qa  is  not  required  in  these 
experiments. 

A  reasonable  measure  of  the  width  of  the  frontal  zone  is  given  by 

d  =  R(-Y,  z,  o)  -  9(Y,  z,  o)   , 

iV^Ux  '''"' 

where  hQ/by   is  approximated  by  a  one-sided  difference.   Figure  1  shows  the 
variation  of  d  for  3  different  experiments  at  the  lowest  level  (z  =  167  m) , 
The  lowest  curve  which  is  taken  from  W72  is  for  the  frictionless  case 
(A  =  0)  with  /ly  =  20  km.   In  this  experiment  no  steady-state  is  reached 
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although  d  is  bounded  by  £y  because  of  the  truncation  error.   In  W72  it 
was  estimated  that  a  discontinuity  would  occur  at  t  =  1.75  days  if  there 
were  no  numerical  truncation  error.   The  middle  curve  corresponds  to  the 
experiment  with  A  =  3  x  lO^m^sec"^  and  Ay  =  20  km,  and  shows  that  the 
solution  reaches  a  steady-state  at  2  days  and  is  close  to  the  limiting 
scale  within  1.5  days.   The  upper  curve  is  for  A  =  5  x  10^  m^sec"^  and 
Ay  =  60  km.   This  experiment  shows  the  same  behavior  as  with  the  smaller 
diffusion  coefficient  except  that  the  limiting  scale  is  much  larger. 

Figure  2  contains  cross  sections  of  0*  =  9  -  0,  ,  where  9  is  the 
horizontally  averaged  initial  potential  temperature.   Only  the  lower 
half  plane  is  shown  because  the  solutions  of  this  section  satisfy  the  sym- 
metry condition  which  is  given  by  Eq.  (5.4)  in  W72.  Also,  the  outer 
portion  of  the  domain  between  Iyj.=  1200  km  and  [y]  =  1800  km  is  not 
shown.   The  solutions  shown  in  the  remainder  of  this  section  are  for 
t  =  4  days.   The  solutions  in  Fig.  2  may  be  compared  with  the  initial 
conditions  and  the  frictional  solutions  which  are  shown  in  Figure  4 
of  W72.   The  solution  for  A  =  3  x  10  m^sec"^  is  very  similar  to  the 
frictionless  solution  at  t  =  1.25  days  except  that  the  frictional  case 
has  a  larger  tilt.   The  solution  for  A  =  5  x  lO^m^sec"''-  is  similar  to  the 
solution  for  A  =  3  x  10^in2  g^^l  except  that  the  gradients  are  much 
weaker  and  the  frontal  zone  has  slightly  less  tilt. 

Figure  3  contains  cross  sections  of  u  for  the  same  cases  and  these 
can  be  compared  with  Figure  5  in  W72.  As  noted  above,  the  solution  for 
A  =  3  X  lO^m^sec"-'-  is  very  similar  to  the  frictionless  solution  at  t  =  1.25 
days,  except  that  the  tilt  is  greater  for  the  frictional  case.   The  solution 
for  A  =  5  X  lO^m^sec"^  is  similar  to  the  solution  A  =  3  x  lO^m^sec"^ 
except  that  the  velocities  are  much  smaller  and  the  gradients  are  much 
weaker . 

Figure  4  summarizes  all  the  numerical  experiments  which  contain 
only  horizontal  diffusion.   The  figure  contains  the  frontal  scale  at 
z  =  167  m  as  a  function  of  the  diffusion  coefficient  on  a  log-log  plot. 
The  dots  represent  experiments  with  Ay  =  60  km  and  the  circles  represent 
experiments  with  Ay  =  20km.   For  comparison,  let  us  derive  the  quasi- 
geos trophic  solution  at  z  =  0.   The  steady-state  first  law  of  thermodynamics 
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at  the  surface  is 


-DyBa'=Aa29'   ^  (4^3) 

^    ~^ 

The  solution  to  this  equation  is       y 

e'Cy,  o)  =  (e'(=o,o)  -  fl'(-oo,o)|(2An/D)^'^'^  expj-pz?.)  dy  .  (4.4) 

o       "*• 

When  this  form  is  substituted  into  the  definition  we  obtain 

d  =  (2Arr/D)^^   ,  (4.5) 

which  is  shown  in  Figure  4.   The  plotted  data  approach  the  line  (4,5) 
as  A  becomes  large,  which  is  expected  because  the  larger  scales  involved 
should  lead  to  more  nearly  quasi-geostrophic  solutions  (Williams,  1967). 
As  A  is  decreased  the  points  fall  fiartiier  below  the  curve  although  the 
slope  is  not  much  greater  than  ^ .   This  is  the  effect  of  the  advections 
of  temperature  and  momentum  by  the  divergent  part  of  the  wind.   These 
effects  are  neglected  in  the  quasi-geostrophic  equations.   For  the  smallest 
values  of  A  the  plotted  data  again  approach  and  cross  the  quasi-geostrophic 
curve.   This  is  a  consequence  of  the  truncation  error  which  arises  when 
the  scale  becomes  very  small. 

5.   Numerical  experiments  with  vertical  diffusion 

In  this  section  we  examine  the  numerical  solutions  which  include  only 

the  vertical  diffusion  so  that  A  =  0,  Both  vertical  diffusion  coefficients 

2-1       2-1 
are  varied  over  a  range  from  1,25  m  sec   to  15  m  sec   .  The  grid  increments 

for  all  of  these  experiments  are  Ay  =  60km  and  Az  =  200m,   This  value  of 
the  vertical  increment  was  chosen  to  resolve  better  the  vertical  structure 
of  the  frontal  zone  and  the  surface  friction  layer.   In  these  experiments 
the  convective  adjustment  process  Q^  is  required  at  times  near  the  surface 
on  the  cold  air  side  of  the  front. 

In  most  of  the  experiments  Cm  ®^d  Cq  are  not  equal  because  they  are 
normally  different  in  the  lower  troposphere  due  to  the  influence  of  thermal 


-10- 


stratification.   Figure  5  shows  the  variation  of  the  frontal  scale  for 
3  combinations  of  C^  and  Cg.   These  curves  may  be  compared  with  Figure  1 
which  gives  solutions  with  horizontal  diffusion  only.   The  solutions  in 
Figure  5  all  show  oscillations  which  were  not  present  in  the  horizontal 
diffusion  cases.   These  oscillations  arise  from  the  development  of  the 
surface  boundary  layer  which  is  not  present  in  the  initial  conditions. 
The  interior  fields  must  adjust  to  the  no-slip  velocity  condition  and  to  the 
insulated  potential  temperature  condition  at  the  lower  boundary.  The 
delay  in  frontogenesis  in  the  first  few  hours  between  Figures  5  and  1 
is  also  due  to  the  adjustment  process.  At  t  =  4  days  the  upper  2  curves 

in  Figure  5  are  nearly  steady  while  the  lowest  curve  (Cm  =  10m2sec"^, 

2   -1 
Cq  =  1.25m  sec  )  still  contains  oscillations.  We  will  compare  the 

various  solutions  in  this  section  at  t  =  4  days  although  those  experiments 
which  have  a  small  Cg  may  not  be  quite  steady. 

Figure  6  contains  cross  sections  of  0   for  the  3  experiments  at 
t  =  4  days.   The  solutions  in  this  section  do  not  have  the  vertical  sjmnnetry 
properties  which  were  noted  in  the  previous  section  because  the  upper  and  lower 
boundary  conditions  are  not  equal.  Nevertheless  we  display  only  the 
solutions  for  the  lower  half  of  the  domain  because  we  are  concerned  mainly 
with  fronts  near  the  surface  and  also  because  the  solutions  near  the  upper 
boundary  are  not  very  realistic.   The  solutions  in  Figure  6  may  be  compared 
with  the  horizontal  diffusion  solutions  which  are  given  in  Figure  2, 
Comparison  of  these  figures  shows  some  important  differences  between 
steady-state  fronts  in  which  only  vertical  diffusion  is  present.   The 
fronts  with  vertical  diffusion  all  show  a  greater  tilt  than  those  with  only 
horizontal  diffusion.   The  solution  with  the  smallest  vertical  heat  dif- 
fusion coefficient  (Cm  =  10  m^sec"-'-,  Cg  =  1,25  m^sec"-"-)  has  a  large  tilt 
in  the  lower  2  kilometers  of  the  atmosphere.   The  vertical  diffusion  in 
these  experiments  causes  the  large  gradient  region  to  extend  higher  into 
the  atmosphere  than  for  those  experiments  with  only  horizontal  diffusion. 
It  is  noted  that  ^fl '/Sz  is  negative  near  the  lower  boundary,  and  this  comes 
from  the  insulated  boundary  condition  at  z  =  0: 

hz  Sz   ^z 

since  Sflj/Bz  is  positive. 
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Figure  7  contains  cross  sections  of  the  u  field  for  the  3  experiments 
at  t  =  4  days.   These  solutions  can  be  compared  with  the  horizontal  dif- 
fusion solutions  which  are  shown  in  Figure  3.  All  of  the  solutions  show 
the  presence  of  the  surface  friction  layer  with  the  no-slip  boundary 
condition.   The  thickness  of  the  surface  boundary  layer  appears  to  vary 
in  proportion  to  Cm^^  ^^  would  be  predicted  by  Ekman  theory.  When  these 
fields  are  compared  with  the  temperature  fields  in  Figure  6,  it  can  be 
seen  that  the  regions  of  large  cyclonic  shear  and  large  vertical  shear 
occur  in  the  zones  of  large  temperature  gradient.  As  would  be  expected, 

the  maximum  speed  occurs  in  the  experiment  with  the  minimum  value  of  the 

2-1  2-1 

momentum  diffusion  coefficient  (Cm  =  1.25  m  sec   ,  Cn  =  5  m  sec  ). 

Figure  8  summarizes  all  the  numerical  experiments  which  contain  only 

vertical  diffusion.   The  data  for  each  experiment  at  t  =  4  days  and  z  =  100  m 

are  plotted  as  a  function  of  C^  and  Cn.   The  number  above  each  point  is 

the  frontal  scaled  and  the  number  below  each  point  is  the  frontal  location 

y^.  A  striking  feature  of  the  diagram  is  the  tendency  for  d  to  remain 

relatively  constant  as  C^   is  increased  when  Cn  is  held  fixed.   Conversely 

d  increases  rapidly  with  Cq  when  C^   is  held  fixed.   The  quantity  y^  is 

a  measure  of  the  mean  frontal  tilt  because  in  the  center  of  the  domain 

(z  =  4.5  km)  the  frontal  zone  is  located  near  y  =  0  (see  Figures  2  and  6). 

As  noted  above  some  experiments  have  a  much  larger  tilt  in  the  lower  2 

kilometers.   Figure  8  shows  a  slow  increase  in  the  tilt  as  Cm  is  increased 

for  a  fixed  C...  When  C^  is  held  fixed  the  tilt  decreases  more  rapidly 

as  Cq  is  increased.   The  experiments  for  the  special  cases  C  =  0  and 

2   -1 

Cq  =  0  are  not  included  in  the  figure.  The  experiment  for  Cjjj  =  10  m  sec 

and  Cp  =  0  contains  large  special  oscillations  while  the  solution  for  Cm  =  0 

2-1 
and  Cq  =  5  m  sec    contains  large  time  fluctuations. 

Some  of  the  gross  features  of  Figure  8  may  be  interpreted  as  follows. 

Consider  the  behavior  for  fixed  C-,  As  Gjj^  is  increased  we  expect  the 

vertical  wind  to  decrease  and  through  the  thermal  wind  equation  this  should 

decrease  the  horizontal  temperature  gradient.   There  is,  however,  an 

opposing  effect  in  the  surface  boundary  layer.   In  the  region  of  cyclonic 

vorticity  we  expect  horizontal  mass  convergence  in  the  boundary  layer; 
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this  should  occur  even  if  the  boundary  layer  is  not  of  the  pure  Ekman  type 
(Mak,  1972).   This  convergence  will  tend  to  build  up  the  temperature 
gradients  within  the  friction  layer.   This  effect  is  opposed  in  part  by 
the  vertical  diffusion  of  heat  which  is  proportioned  to  Cq.  As  Cjjj  is 
increased  the  boundary  layer  will  become  thicker,  as  was  pointed  out  earlier 
in  this  section.   This  increase  in  the  thickness  of  the  friction  layer 
will  tend  to  increase  the  temperature  gradients  as  long  as  Cq  is  held 

u 

fixed.   The  numerical  solutions  show  that  this  frictional  convergence 
effect  on  the  temperature  gradient  is  approximately  balanced  by  the  momen- 
tum diffusion  effect  through  the  thermal  wind.  Therefore  d  remains 
nearly  constant  as  C^   is  increased  with  Cq  held  fixed. 

When  C  is  held  fixed  and  Co  is  decreased  the  vertical  scale  of  the 

m  0 

frontal  zone  becomes  very  small.   This  can  be  seen  in  Figure  6  if  the  lower 
2  cross  sections  are  compared.   The  horizontal  temperature  gradient  is 
somewhat  constrained  through  the  thermal  wind  relation  by  vertical  momentum 
diffusions.   The  decrease  in  the  vertical  scale,  which  comes  from  a  de- 
crease in  Cq,  must  be  associated  with  an  increase  in  the  tilt  of  the 
frontal  zone  if  the  change  in  the  horizontal  gradient  is  to  be  minimized. 
This  behavior  is  clearly  present  in  Figure  8. 

6.   Numerical  experiments  with  both  horizontal  and  vertical  diffusion 

In  this  section  we  will  examine  the  solutions  which  contain  both 
horizontal  and  vertical  diffusion.   Individual  solutions  will  not  be 
shown  because  they  are  essentially  combinations  of  the  solutions  shown  in 
earlier  sections.  The  principal  features  of  the  solutions  are  summarized 
in  Tables  1  and  2.  All  solutions  have  space  increments  of  ^  =  60  km 
and  t:^  =   200  m.   The  data  are  extracted  at  t  =  4  days,  and  all  of  the 
data  in  the  tables  are  for  the  lowest  level  (z  =  100  m) . 

Table  1  contains  the  frontal  scale  d  for  various  combinations  of  the  dif- 
fusion coefficients.   The  first  column  gives  the  values  for  aero  vertical 
diffusion  and  the  first  row  gives  the  values  for  zero  horizontal  diffusion. 
The  smallest  value  of  d  occurs  at  the  lowest  level  for  each  experiment 
except  for  the  two  with  footnotes.   In  the  two  exceptional  cases,  the 
minimum  occurs  at  the  second  point  and  this  is  caused  by  the  vertical 
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diffusion  of  momentum.   It  may  be  seen  that  the  value  of  d  for  each 
experiment  is  a  little  less  than  the  sum  of  the  values  obtained  from  the 
experiments  where  one  of  the  coefficients  is  zero. 

Table  2  contains  the  frontal  locations  y^  for  various  combinations  of 
the  diffusion  coefficients.   In  the  first  column,  which  is  for  zero  ver- 
tical diffusion,  the  tilt  decreases  slightly  with  increasing  horizontal 
diffusion  coefficient.   The  combined  experiments  all  show  values  which 
are  closer  to  the  A  =  0  than  to  the  Cjq  =  Cg  =  0  solutions.  However, 
for  each  combination  of  C  and  Cq,  the  tilt  does  decrease  as  A  is  in- 
creased. We  conclude  that  the  frontal  scale  is  determined  by  both  dif- 
fusion processes,  but  that  the  frontal  tilt  is  determined  principally 
by  the  vertical  diffusion, 

7.   Conclusions 

The  numerical  model  of  Williams  (1972)  is  modified  in  this  study  to 
include  horizontal  and  vertical  diffusions  of  momentum  and  heat.   The 
turbulent  diffusions  are  represented  with  constant  coefficients.  An  Ekman 
layer  is  added  to  the  basic  deformation  field  and  some  other  teirms  which 
are  related  to  the  modified  deformation  field.   The  numerical  solutions 
become  quasi- steady  within  1-2  days.   These  results  suggest  that 
atmospheric  fronts  form  in  a  fairly  short  period  of  time  (1-2  days) 
and  then  move  in  a  state  of  quasi-balance  where  the  frontogenetic  processes 
are  balanced  by  diffusion  processes.   The  time  scale  of  this  second  phase 
would  be  related  to  the  time  scale  of  the  deformation  field  which  caused 
the  frontogenesis. 

The  solutions  with  only  horizontal  diffusion  are  similar  to  the  transient 
solutions  obtained  in  W72  at  particular  times.  The  horizontal  scale  at 
the  lowest  level  follows  approximately  an  A^  dependence.   The  values  obtained 
are  less  than  the  exact  quasi-geostrophic  solutions  at  z  =  0  and  they  approach 
the  latter  as  A  increases.   For  reasonable  values  of  A  the  scale  is  close 
to  the  quasi-geostrophic  value,  but  the  structure  is  very  different  since 
the  quasi-geostrophic  solutions  contain  no  vertical  tilt. 

The  solutions  with  only  vertical  diffusion  are  strongly  dependent 
on  the  ratio  of  C^  to  Cg,  All  of  these  solutions  show  a  deeper  region 
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of  large  gradient  than  for  the  solutions  with  only  horizontal  diffusion. 
When  Cq  is  held  fixed,  the  scale  d  shows  little  change  as  C^   is  increased. 
On  the  other  hand,  when  Cm  is  held  fixed,  d  increases  rapidly  with  increasing 
Cq.   The  frontal  tilt  increases  when  Cm  is  held  fixed  and  Cq  is  decreased. 
The  solutions  for  C^j^  =  0  or  Cq  =  0  are  not  acceptable  steady  solutions. 

The  solutions  which  include  both  horizontal  and  vertical  diffusion 
are  combinations  of  the  separate  solutions.   The  values  for  d  are  a  little 
less  than  the  sum  of  the  values  from  each  of  the  diffusion  processes. 
The  tilt  of  the  fronts  is  principally  determined  by  the  vertical  diffusion 
effects. 

The  steady-state  frontal  solutions  obtained  in  this  paper  are  reasonable 
with  respect  to  observed  fronts.   The  vertical  diffusion  terms  are  clearely 
required  to  give  the  proper  structure  near  the  surface.  Both  Cm  and  Co 

must  be  included  although  a  value  C  that  is  somewhat  less  than  Cm  n»ay  be 

B 

appropriate.   The  importance  of  the  horizontal  diffusion  processes  can 
probably  be  determined  only  from  turbulence  studies  in  frontal  zones. 

The  solutions  found  in  this  paper  would  be  affected  if  nonlinear 
diffusion  coefficients  were  used,  but  we  expect  the  same  general  type  of 
solution.  However,  if  C  were  allowed  to  decrease  above  the  friction  layer, 
the  maximum  gradient  might  occur  higher  in  the  atmosphere. 

The  tilt  of  the  solutions  obtained  would  be  modified  if  the  whole 
system  were  allowed  to  have  a  net  translation.   Then  it  would  be  possible 
to  distinguish  between  cold  and  warm  fronts.   The  addition  of  latent 
heat  should  be  important  in  determining  the  structure  of  many  fronts. 
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Table  1.   Values  of  d  in  kilometers 


A    ^\ 

0,0 

10,  1.25 

10,  10 

1.25,  5 

0 

172 

399 

241 

3  X  10^ 

94 

239 

450 

291 

10^ 

181 

321 

523 

370 

5  X  10^ 

492 

666^ 

2 
828 

644 

1  Minimum  at  z  =  300  m  of  631. 

2  Minimum  at  z  =  300  m  of  827, 

Table  2.   Values  of  y^  in  kilometers 


A  ^\ 

0,0 

10,  1.25 

10,  10 

1.25,  5 

0 

—  M  — «• 

-1010 

-820 

-720 

4 
3  X  10 

-545 

-995 

-805 

-715 

10^ 

-540 

-985 

-755 

-700 

5  X  10^ 

-510 

-985 

-715 

-640 
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FIGURE   CAPTIONS 

Figure  1.   Time  variation  of  the  frontal  scale  at  z  =  167  m  for  experiments  with 
only  horizontal  diffusion. 

—  5   2   -1 

Figure  2,   Cross  sections  of  9'  =  9  -  Bj  for  A  =  5  x  10  m  sec   and 

A  =  3  X  10^  m^sec"!  at  t  =  4  days. 

5  2-1  4  2-1 

Figure  3.   Cross  section  of  u  for  A  =  5  x  10  m  sec   and  A  =  3  x  10  m  sec 

at  t  =  3  days. 

Figure  4.   The  variation  of  the  frontal  scale  at  z  =  167  m  as  a  function  of  the 
horizontal  diffusion  coefficient. 

Figure  5.   Time  variation  of  the  frontal  scale  at  z  =  100  m  for  experiments 
with  only  vertical  diffusion. 

Figure  6.   Cross  section  of  9'  for  3  combinations  of  C  and  C^  at  t  =  4  days. 

m      9 

Figure  7.   Cross  sections  of  u  for  3  combinations  of  C  and  C^  at  t  =  4  days. 

Figure  8.   The  frontal  scale  and  frontal  location  as  functions  of  C  and 
C  .   These  data  are  for  z  =  100  m  and  t  =  4  days. 
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Fig.    6 
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Fig.    7 
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Fig.    8 
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